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Abstract. O(N) invariant vector models have been shown to possess non-trivial scaling 
large iV limits, at least perturbatively within the loop expansion, a property they share 
with matrix models of 2D quantum gravity. In contrast with matrix models, however, 
vector models can be solved in arbitrary dimensions. We present here the analysis of 
field theory vector models in d dimensions and discuss the nature and form of the critical 
behaviour. The double scaling limit corresponds for d > 1 to a situation where a bound 
state of the iV-component fundamental vector field <p, associated with the <p 2 composite 
operator, becomes massless, while the field <p itself remains massive. The limiting model 
can be described by an effective local interaction for the corresponding O(N) invariant 
field. It has a physical interpretation as describing the statistical properties of a class of 
branched polymers. 

It is hoped that the O(N) vector models, which can be investigated in their most 
general form, can serve as a test ground for new ideas about the behaviour of 2D quantum 
gravity coupled with d > 1 matter. 



1. Introduction 

Vector models have, in zero dimension, in the large N limit, critical points and scaling 
behaviours [|],|2]||] that are reminiscent of those of matrix models relevant to 2D quantum 
gravity While matrix models are associated with surfaces because the large iV 

limit selects planar Feynman diagrams, vector models are associated with one-dimensional 
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chains because only so-called "bubble" diagrams survive for iV large. The correspond- 
ing critical points seem indeed to have a natural physical interpretation in the theory 
of branched polymers or filamentary surfaces 0,0. We report here the extension of the 
analysis to vector models in quantum mechanics and d dimensional field theories [jSJ. 

These peculiar critical points of vector models have the following origin: In the large N 
limit, vector models are reduced to perturbation theory around minima of an effective po- 
tential resulting from the competition between an attractive potential and a measure term 
which suppresses contributions of small fields. A singularity occurs when the minimum of 
the effective potential disappears. In matrix models similarly there is a competition be- 
tween the potential and a measure term which produces a repulsion between eigenvalues, 
although the detailed mechanism is slightly different. 

Vector and matrix models share also other properties: Perturbation theory is always 
divergent @ and half of the vector models (like for matrix models) are ill-defined beyond 
perturbation theory because perturbation theory is non-Borel summable and the corre- 
sponding potential is not bounded from below. This is a serious problem because in both 
cases this class includes the models with positive weights, the only ones which have a 
simple physical interpretation. 



2. O(N) -invariant integrals 

Branched polymers and O(N) -invariant models. Consider the integral 

Z = J d^xda exp (-x 2 - ax 2 - Na 2 /\) . 

We can integrate over the iV-component vector x: 

Z = J da exp [-N (a 2 /A + \ ln(l + a))] . 

If we then expand the partition function in powers of A, we observe that the coefficients 
of the expansion are sums of Feynman diagrams which can be interpreted as some kind 
of branched polymers. These polymers are weighted in particular by a factor X 1 where 
/ counts the number of links (generated by the a-propagator) and a factor N X ~ L where 
L is the number of loops. As in matrix models the continuum limit is reached when A 
approaches the singularity of Z closest to the origin, emphasizing the long chains. In the 
large N limit only tree-like chains are selected. We shall verify that there exists a double 
scaling limit in which chains with an arbitrary number of loops contribute. 
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To see the connection with 0(iV)-invariant vector models we instead integrate over a 
and obtain 

Z = Jd N x exp[-x 2 + A(x 2 ) 2 /4] • 

General O(N) -invariant vector models. We immediately consider more general O(N)- 
invariant vector models with a partition function given by the iV-dimensional integral: 

d ff xe-^, (2.1) 

where V(a) a polynomial. Integrating over angular variables we remain with an integral 
over a = x 2 

ZoC r^ e -mV(a)-^na]^ (2 _ 2) 

Jo °~ 

With this more general integrand we can vary independently the weights associated with 
the different vertices. The large N limit can be evaluated by steepest descent. The saddle 
point a s is given by: 

2V'(a s )a s = 1. (2.3) 

Critical points are points where the derivative of V'{a)a at a s vanishes so that the result, 
at leading order, is no longer given by the gaussian approximation. Let us assume more 
generally that p — 1 successive derivatives vanish. The critical potential V = V c then 
satisfies 

2V^(a)a-l = 0((a s -a) p ), V c (a) - \ In a = const. + 0((a s - a) p+1 ). 

We see that the values of a — a s contributing to the integral are of order jVVfr+i). Let 
us add to the critical function V c the set of relevant perturbations (in the sense of critical 
phenomena) 

p-i 

V(o-) = V c (o-) + Y,v q (°s-o-) q , 

q=l 

(the term q = p can be eliminated by a shift of a). Rescaling a s — a into (a s — a)N~ x / ( p+1 ) 
we see that the scaling region corresponds to take v q = 0(N < ^ q ~ p ~ 1 ^^ p+1 ^). Setting then 

a s - a = zN- 1 /^, u q = N^^l^v, , 

we find the universal scaling free energy F = \nZ: 

F{u q ) ~ In J dz exp (^-z p+1 + ^u g z 9 j • (2.4) 
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Note that the case p odd corresponds to convergent integrals while the case p even can be 
only reached by analytic continuation and yields complex functions. This is in particular 
the case for the p = 2 critical point which is the only one which has a simple statistical 
interpretation because all weights are positive. A similar situation is found in matrix 
models, where in particular the unitary models are all unstable. 



3. O(N) invariant potentials in quantum mechanics 

We consider now the partition function Z 



[dx(f )] e- 



-S(x) 



(3-1) 



for the O(N) invariant action 



S — N dt 



i(x) 2 + ix 2 + ^) : 



(3.2) 



At the end of the section we shall briefly examine the problem of the existence of a large 
N scaling limit, using radial coordinates, from the Schrodinger equation point of view 
(the only method available in the matrix case). However, because we want to extend the 
analysis to field theory, we use here a more general method, which directly shows the 
relation between the vector model and a model of branched polymers embedded in one 
dimension. We start from the identity 



exp 



NX 



dt (x 2 ) : 



oc J [dcr( £)] exp 



(3.3) 



Introducing identity ( |3.3|) into the path integral, we can integrate over x and find 



Zoc J [da(t)]e- Scff(fj) , 



with 



S e s(cr) 



N 



dt 



° 2 {t) 
2A 



+ tr In (-d 2 + 1 + a) 



(3.4) 



(3.5) 



The dependence of the partition function on iV is now explicit and the partition function 
can be interpreted in terms of a branched polymer model, the a internal lines representing 
the polymer links. 



In the large N limit the path integral can be calculated by steepest descent, 1/N 
playing the role of % in the classical limit. We look for a saddle point a(t) = a constant 
which is then a solution of the equation: 



a 
A 



1 

2^ 



a 
A 



1 







u 2 + l + a A 2tJT+~g 
A critical point is found when two different solutions meet. The condition is 



(3.6) 



d 

d7 



a 
A 



1 



l 



+ 



l 



A 4(1 + (T) 3 / 2 



0, 



(3.7) 



This equation implies that A is negative, and thus the potential is not bounded from below. 
The analysis which follows has thus no meaning beyond perturbation theory. Although 
the sum of the perturbative expansion can be defined by analytic continuation, the result 
is then complex and thus physically not acceptable. 
With our normalizations the explicit values are 



A c = -4/3 3/2 < , a c = -2/3. 



(3.8) 



The 1/N -expansion. To generate the 1/N expansion we have first to calculate the 
ex-propagator. It is convenient to set: 

= 1 + a fi 2 c = 1/3. 

The Fourier transform A a (u) of the cx-propagator is then: 



1 1 

- — h — — 



du/ 



A 2tt J ( w /2 + /i 2) hu-u'f + n 



(3.9) 



The criticality condition implies that the inverse propagator vanishes at oj = 0. Note that 
since fi c does not vanish the inverse-propagator remains a regular function of cu 2 , which 
behaves like uj 2 for uj small. In higher dimensions in the critical limit a would become a 
massless field . In one dimension, however, quantum fluctuations always lift the degeneracy 
of the classical ground state. To understand the critical limit we perform a local (small u 
in Fourier transform) and small a — a c expansion of the action. We then rescale time and 
a: 

t i-> tN 1/5 <=^> u h-> uN~ 1/5 , a-a c ^&N~ 2/5 , (3.10) 
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to render the coefficients of a 3 and (a) 2 in the action N independent. We note that terms 
with higher derivatives or higher powers of a are then suppressed for iV large. Therefore 
at leading order for iV large (and after some additional finite rescaling) the action takes 
the form: 



S(a) 



dt 



2 \~dtj + T 



(3.11) 



The action is not bounded from below and the path integral can be defined only by analytic 
continuation. The corresponding the hamiltonian H is: 

H = N~ 1 / 5 (-I(d/da) 2 + |a 3 ). 



The scaling region. For A close to A c the most relevant new interaction term is the 
term linear in a. After the rescaling ( |3.1U ) it gives an additional contribution to the action 
proportional to (A — A c )iV 4 / 5 cr. The scaling region is thus defined by letting iV go to infinity 
and A — A c to zero at 

(A - A c )iV 4 / 5 



u 



fixed. The hamiltonian relevant to the scaling limit is then: 



H = iV~ 1//5 



1/dV ^ 

— — +ua-\ 

2 1 da J 3 



(3.12) 



We have found a scaling regime analogous to the one observed in d < 1 matrix models. In 
the d = 1 matrix model, instead, the situation is more complicated |l(| because logarithmic 
deviations from a simple scaling law are found, a situation we shall meet in the d = 2 vector 
model. 

Generalization. The previous analysis can be extended to more general (9(iV)-invariant 
potentials. Let us consider the action: 



S = N dt 



i(x) 2 + F(x 2 ) 



(3.13) 



The large iV expansion is generated by standard techniques (see for example |TT|). We 
introduce a Lagrange multiplier p to impose the constraint a = x 2 . The action then takes 
the form 



S(x, p,a) = N dt \ (x) 2 + V{a) + \p (x 2 - a) 



(3.14) 



We integrate over x to obtain: 

S(p, a)=N^J dt [V(a) - \pa\ + ±tr In (-d 2 + p) j . (3.15) 

In the large N limit the path integral can again be calculated by steepest descent. We 
look for two constants p, a solutions of 

// x 1 i 1 f du i 1 , 

V'(a) -hp = 0, ~ h<r + — / -5 = -h<?+ -r^ = • 3 - 16 

v ' 2 2 4tt J u 2 + p 2 4^/p v ; 

When the determinant of partial derivatives of this system vanishes we find a critical point: 

detS^detf^ _ri/ 2 )=0, 
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and thus 

8a 2 V'(a) = 1, ^V^a) = -1. 

The matrix 5 , ^ 2 - ) is the inverse propagator at oj = 0. These conditions thus imply that 
the propagator has a pole at u> = 0. However only the linear combination a of p and a 
which corresponds to the eigenvector of with zero eigenvalue is singular. The non-zero 
mode can be integrated out and we find an effective action for a. The main difference with 
the case of the quartic interaction is that we can now adjust the original potential V in 
such a way that the coefficient of the leading cr 3 interaction vanishes, or more generally 
all interactions up to the a v vanish. We then perform a local expansion of the action and 
rescale time and a to render the coefficients of (a) 2 and a p+1 N independent: 

t tjV(P-l)/(P+3) ^ u _> UN -( P -D/(P+S) ^ ~ _ ~ c _> ~ Ar -2/( P+ 3)_ (3 _ 17) 

It is easy to verify that interactions containing derivatives or higher powers of <r are then 
suppressed for large. 

To describe the scaling region we can add to the critical potential a set of relevant 
terms characterized by parameters v q . The scaling limit is then obtained by keeping the 
products u q , 



u 



q = N (W-W(r+*) Vqi q =l,..., p -l, 



fixed and the corresponding scaling hamiltonian is: 



H = A^" (p_1)/(p+3) 



l/d\ 2 ^ 

2UJ + ^TT + g^7 



(3.18) 



Again only half of the models corresponding to p odd are stable. 

Hamiltonian formalism. For A" large the zero angular momentum hamiltonian Hq can 
be written (after factorizing r^ -1 )/ 2 in the wave function): 

i fa-* 



where r = |x| and for iV large W(r) is related to the potential V(r) by 

W(r) = ± + V(r). 
In the case of the anharmonic oscillator (|3.2|) for example 
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For A" large, the eigenvalues can be calculated by expanding perturbation theory around 
the classical minimum r c of the potential W. A critical potential is defined by the condition 
that the second derivative of W also vanishes: 

W'(r c ) = W"(r c ) = 0. (3.20) 

If moreover p derivatives of W vanish the leading order hamiltonian is 

H c = NW(r c ) + (^^VcXr - r c )^\ 

which leads after rescaling to a critical contribution to the ground state energy of order 
n -(p-i)/(p-3) in agreement with (|3TT8|) . 



In the special case Q conditions flOOp yield r 2 = 3 1 / 2 /2, A = -4/3 3 / 2 . We 
recognize the critical value ( |3.8|) of A. 



4. Field theory: Critical points in the large limit 

We come now to the most interesting case, field theory in d > 1 dimensions, which we 
expect, according to the previous analysis, to correspond to branched polymers embedded 
in d- dimensional space. 

In contrast to quantum mechanics (d = 1) we expect now the a-field to remain massless 
even after taking into account the successive corrections of the large expansion. The 



8 



field a, which is equivalent to the composite (ft 2 field, is associated with a massless bound 
state of the field (ft, which itself remains non-critical. 

We discuss the problem only in the special case of the (4> 2 ) 2 interaction, the extension 
to more general cases being straightforward. We consider the partition function: 



[d<^(cc)] e" 



-S(4>) 



where S((ft) is the O(N) symmetric action: 

S(<ft)=N J {l[d,<P(x)} 2 + 1 -r<p 2 (x) + ±[ ( p 2 (x)} 2 }a d 



(4.1) 



(4.2) 



A cut-off of order 1, consistent with the symmetry, is now implied. For example we can 
assume that the inverse propagator has higher order derivative terms and we have explicitly 
written in action (\i.2\j only the two first terms in a local expansion (in Fourier space small 
momentum expansion). In particular, the parameter r is then the value of the inverse 
propagator at zero momentum. 

We use the same algebraic identity as in the quantum mechanical case 

NX 



exp 



d a x 



oc J [da(x)} exp 



N 



d x — 6 Z 

V4A T 



Introducing this identity into ([4JJ) and integrating over (ft we find 

Z = j /W)]e-^> 



with: 



N 
~2 



f a 2 (x) 

/ — K - 1 d d x + tr ln(-A + r + a(x)) 
J 2A 



(4.3) 



(4.4) 



(4.5) 



expression which shows that the vector model is now related to a branched polymer problem 
in d dimensions. 

The large N limit. In the large iV limit the functional integral can be calculated by 
steepest descent. We look for a uniform saddle point a(x) = a. The saddle point equation 
is: 

it i f d v 

0. (4.6) 



a i 

A + (2^p 



p z + r + a 



Let us introduce the parameter m: m 2 = r + a, which is, in the large N limit, the mass of 
field (ft. Eq. (|4.6|) can then be written: 



(m 2 — r) 



A 



1 



d d p 
p 2 + m 2 



0. 



(4.7) 



The solution is singular when the derivative with respect to a or m 2 vanishes. This yields 
the equation for the critical point: 

— + -r / — ~ = . 4.8) 

Eqs. (|4.6| ), (|4.8| ) define, at r fixed, critical values of A and a. The critical value A c of A 
is again negative. Criticality eq. ( }4.8[ ) implies that the a-propagator has a pole at zero 
momentum. The a-field becomes massless while the </>-field remains massive. 



5. The double scaling limit 

We have studied the large N limit. We now look for a scaling limit. Since the a field 
is a one-component field it can remain critical for d > 1 even in presence of interactions. 
Still perturbation theory is IR divergent for dimensions d < 6. We therefore add a relevant 
perturbation proportional to A — A c which provides the theory with an IR cut-off. We then 
look for the most IR divergent terms in perturbation theory: This is a problem standard 
in the theory of critical phenomena |[11||, the bare mass squared which is proportional to 



VA C — A playing the role of a deviation from the critical temperature. The effective action 
for the a-field is non-local and contains arbitrary powers of the field. However, because 
the (p-field is not critical we can again make a local expansion. Standard arguments of 
the theory of critical phenomena tell us that the most IR divergent terms come from 
interactions without derivatives and with the lowest power of the field. Here the leading 
interaction is proportional to a 3 . To characterize the IR divergences of the perturbative 
expansion in powers in 1/N we again rescale distances and field a — a c : 

a - a c oc aA (2 " d)/2 iV- 1/2 , x i-> Ax, (5.1) 

where A will play the role of a cut-off. The effective action at leading order, after some 
additional finite renormalizations, is 



S eff (cr) = J d d x 



A(6-d)/2 



where v oc A — A c and a cut-off A is now implied. 

We first consider dimensions d < 6. The a 3 field theory is then super-renormalizable 
and we fix the coefficient of ex 3 : 

A (6- d )/2 /v /^ = ^ 
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Therefore the cut-off grows with N like jV 1 /( 6_d ). However, unlike the case of quantum 
mechanics, we cannot also keep the quantity vA^ d+2 ^ 2 y/N fixed because the field theory 
has UV divergences when the cut-off A becomes large. 

Dimensions d < 4. It is convenient to examine first dimensions d < 4 because then 
only the field average is divergent. We have to introduce a counterterm which renders 
< &(x) > finite. The renormalized action is 

Seft(a r ) = J d d x [± (d^a r ) 2 + \n 2 ~a 2 r + ^a s r - Cl (A)a r 

in which ay is the renormalized field and \x is a renormalized mass parameter. 

To recover the original action we must eliminate the term quadratic in the field and 
thus shift a r by a quantity a: 

a = -fi 2 /2g 3 . 

Identifying then the coefficients of the linear term we find: 

v A( d + 2 )/ 2 ^N = - Cl (A)-^-. 
For d = 2 only the one-loop diagram is divergent and we obtain 

d(A) = f?-ln(A/A0. 

2,71 

Therefore to obtain a non-trivial scaling limit we have to choose: 



1 

' ~ ~N 



(5.2) 



Note that only the In N/N term is universal, the 1/N term is regularization and renormal- 
ization dependent. 

For d = 3 the two- loop diagram is also divergent. Still for A large the leading contri- 
bution is still given by the one-loop diagram, thus ci(A) oc A and v = 0(1/N). We have 
thus shown that for d < 4 a scaling limit exists which leads to a renormalized a 3 field 
theory. The relation between A — A c and N, however, has itself no longer a simple power 
law form. 

Note finally the similarity between the results for the vector model at d = 2 and the 
matrix model at d = 1 [|10] . 
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Higher dimensions. For 4 < d < 6 the situation is slightly more complicated because 
two counterterms are required, renormalizing < a > and the coefficient of a 2 . The quantity 
v oc (A — A c ) becomes a even more complicated function of N. However, at leading order 
for N large, v still behaves as 1/iV, while naive scaling would have predicted a scaling 
variable aM 6 "^ 4 . 

For d > 6 IR divergences are no longer strong enough to compensate the 1/N factors 
and thus no non-trivial scaling limit can be defined. 

More General Interactions. The method applicable to more general interactions has 
already been explained in the case of quantum mechanics. Because only one mode is 
critical the main effect is to introduce additional parameters in the effective interaction of 
the critical field in such a way that the most IR divergent interactions can be cancelled. 
In the language of critical phenomena we reach multicritical points. In dimensions d < 
2(p + l)/(p — 1) we then generate the renormalized a p+1 interactions provided we again 
choose the parameters of the potential, as functions of N, such that they cancel the UV 
divergences of perturbation theory. 

6. Conclusion, prospects 

We have shown on a few simple examples that vector models have, at least in low 
dimensions, non-trivial large N scaling limits. Obviously more general models can be 
analyzed by similar methods, for example by introducing several vector fields additional 
degrees of freedom for the polymers are generated. Of particular interest is the general- 
ization to models containing fermions (in particular supersymmetric models). One would 
like to find out whether the presence of fermions stabilizes some of the models which are 
unstable otherwise, a very serious problem for the corresponding matrix models of 2D 
gravity. 
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